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Abstract. We use the PHOENIX code package to model the X-ray spectrum of Nova V4743 Sagittarii observed with 
the LETGS onboard the Chandra satellite on March 2003. Our atmosphere models are ID spherical, expanding, 
line blanketed, and in full NLTE. To analyze nova atmospheres and related systems with an underlying nuclear 
burning envelope at X-ray wavelengths, it was necessary to update the code with new microphysics, as discussed 
in this paper. We demonstrate that the X-ray emission is dominated by thermal bremsstrahlung and that the 
hard X-rays are dominated by Fe and N absorption. The best fit to the observation is provided at a temperature 
of T cff = 5.8 x 10 5 K, with L bo i = 50 000 L . The models are calculated for solar abundances. It is shown that 
the models can be used to determine abundances in the nova ejecta. 
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1. Introduction 

Classical novae (CNe) are the third most violent of the 
stellar explosions that can occur in a galaxy (after 7- 
ray bursts and supernovae). About 5 to 10 Galactic 
CNe are discovered per year in the Solar neighborhood. 
CNe participate in the cycle of Galactic chemical evo- 
lution, in which they eject grains and metal-enriched 
gas as a source of heavy elements for the interstellar 
medium (supplementin g those of supernov ae, AGB stars, 
PNe, and WR stars: iGehrz et all Il99s|) . In addition, 
CNe are related to the Super Soft Binary X-ray Sources 
(SSS) which are probably progenit ors of SN la explosions 
jKahabka fe van den Heuvellll997l) and which are thought 
to be the result of thermonuclear explosions in the cores 
of mass-accreting white dwarfs in close binary systems 
(Hillebrandt fe Niemeve J200(j lLeibundgutl200 1]) . The as- 
sociation of CNe and SSS is important both because of 
the cosmological implications of SN la (currently the best 
standard candles for determining cosmological distances) 
and because SN la are thought to be responsible for the 
abundance of iron in the solar system. 

CN explosions occur in a cataclysmic variable (CV) bi- 
nary star system in which a Roche lobe filling secondary 
supplies hydrogen-rich material that is accreted onto a 



white dwarf (WD) primary. Theoretical studies show that 
the accreted layer grows until it is sufficiently hot and 
dense to initiate a thermonuclear runaway (TNR). The 
evolution of the TNR depends upon the mass and lumi- 
nosity of the WD, the rate of mass accretion, the compo- 
sition of the accreted materi al, and the chemical compo- 
sition in the reacting layers l|Starrfieldlll989T: IGehrz et all 
1998j). Observations imply and theory demands that core 
material be dredged-up into the accreted material and 
the ga ses be thoroughly mixed before being e jected into 
space l|Starrfield et al.lll998llGehrz et aljfl99^ . Thus, the 
chemical composition of the ejected material reflects a 
TNR processed mixture of WD core plus accreted mate- 
rial. Abundance studies of the ejecta carried out with the 
IUE satellite have established that the underlying WDs 
are either carbon-oxyge n (CO) or oxygen, neon, and mag- 
nesiu m (ONeMg) WDs llStarrfield et alll998llOehrz et alJ 
Il998^ . CNe are expected to be the major source of 15 N and 
17 O in the Galaxy and to contribute to the abundances of 
other isotopes in this atomic mass range. 

Previous X-ray studies of CNe in outburst show that 
no other wavelength region provides unambiguous infor- 
mation on the evolution and characteristics of the under- 
lying WD [Drake et al.ll2003t iKrautter et ai1il99rl l2002t 
iNess et al.ll2003j) . For example, the ROSAT studies dis- 
covered a "soft" component in CNe and placed them at, 
or near, the bright end of the class of SSS. The "soft" 
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component has been identified as emission from the hot 
nuclear burning photosphere of the WD and theoreti- 
cal studies show that a determination of the tempera- 
ture and lifetime is a measure of the mass of the WD. 
Stellar atmosphere analyses of the photosphere can pro- 
vide abundances of the material remaining on the WD. 
CNe at optical maximum are luminous, L> 1O 5 L0 and 
evolve to X-ray maxi mum where T B ff ranges from 3 x 10 5 K 
to 6xl0 5 K or higher farautter etal1ll996t ISalman et ail 
ll998HOrio et alJ^nO^lTonsh . 

X-ray emission during the CN outburst can be divided 
into three phases l(Krautterl l2002). First is the early fire- 
ball phase, in which the hot atmosphere of the white dwarf 
is expanding and cooling adiabatically. The peak temper- 
ature of the white dwarf is predicted to e xceed 10 6 K (de- 
pend ing on the mass of the white dwarf: IStarrfield et al. 
Il99(il) so, if detected in this phase, it would be an ex- 
tremely hot source. Nevertheless, the expanding shell cools 
rapidly and becomes more opaque in X-rays in a few 
hours. Therefore, for the first few hours of the outburst, 
after the shell has cooled sufficiently, the nova is probably 
undetectable by an X-ray satellite. However, after a few 
days, the expanding shell starts to become ionized and 
the copious X-rays from the inner layers partially pen- 
etrate the shell and X-ray observations of novae during 
this phase show emission lines from the ex panding gas 
llMukai fc Ishid^nOltlOrio et al.ll2002ll20n3h . This is the 
second phase. 

The third is the constant bolometric luminosity phase 
when the system again becomes observable as a "soft" 
source. This phase occurs because both observational and 
theoretical studies of CNe show that not the complete 
envelope is ejected during the explosive phase of the out- 
burst. Some fraction remains on the WD and rapidly re- 
turns to hydrostatic equilibrium. This material provides 
sufficient fuel so that the enlar ged WD remains hot and 
luminous for months to years i Krautter et al.l Fl996^ . It 
is predicted that the duration of this phase is an inve rse 
function of the mass of the WD ( Krautter et al.lll99^ so 
that a determination of how long a nova is "on" in X-rays 
can provide an estimate of the WD mass. Moreover, since 
the burning WD is hot, with nuclear burning ongoing in 
its surface layers, its characteristics closely resemble those 
of the SSS. 

In this paper we report on the analysis of an observa- 
tion of nova V4743 Sgr with the High Resolution Camera 
and the Low Energy Transmission Grating Spectrograph 
(HRC-S + LET GS) on the CHANDRA X-ray satellite 
(|Ness et al.lr2fin3h . Because the structure of the rekindled 
WD in V4743 Sgr consists of a nuclear burning envelope 
with an very hot stellar atmosphere on top, we have up- 
dated the microphysics and numerics in PHOENIX to ana- 
lyze nova atmospheres at X-ray wavelengths. We describe 
these improvements in section l?"2l PHOENIX is well suited 
to analyze V4743 Sgr, or any SSS, because it has been 
developed specifically to treat radiation transport in an 
expanding, optically thick atmosphere at any tempera- 
ture. We show that we can now treat the formation and 



time evolution of the X-ray radiation in a CN in outburst 
and that it is possible to determine the chemical composi- 
tion of the CN ejecta. We describe the X-ray observations 
of V4743 Sgr in section [21 and end with a Summary and 
Discussion (sections 0] and [5J- 

2. Modeling novae in X-rays 

The nova models calculated with the PHOENIX code are ID 
spherical, expanding, line blanketed, NLTE model atmo- 
spheres. The spherical, co-moving frame radiative transfer 
equation (SSRT) for expanding media is solved for lines 
and continua, cou pled with the NLTE statis tical equi- 
librium equations ijHauschildt fc Baronl 11999^1 . The nu- 
merical solution of the radiation transport and multi- 
level NLTE problems are based on an operator splitting 
(A LI) method as described in d etail bv lHauschildtl i|l992f) 
and lHauschildt fc Baronl (^999). Nova atmospheres are far 
from LTE, therefore, all models have to include NLTE for 
as many atoms and energy levels as feasible. 

The radiative transfer problem is coupled to the en- 
ergy conservation through the equation of radiative equi- 
librium (in the co-moving fra me) and, therefore, a te m- 
perature correction procedure l Hauschildt et alJl2003l) is 
used to iteratively correct the structure so that energy 
is conserved. Mechanical energy sources, due to the ex- 
pansion of the nova envelope, as well as convection are 
negligible for the conditions found in nova atmospheres. 

The density profile and the velocity field of the ex- 
panding medium a re taken from h ydrodynamic simula- 
tions l|Sharalll989t IStarrfieldll989|) . In this case the ra- 
diation transport problem effectively decouples from the 
hydrodynamic equations and the problem is dramati- 
cally simplified. Earlier nova models in other wavelength 
ranges had shown that the re are satisfactory results by th e 
use of these simulations l|Hauschildt fc Starrfieldlll99nl) . 
Therefore we use the results of the simulations as a first 
assumption in this stage of our work. 

The nova atmosphere is approximated by an expand- 
ing but stationary (in time) structure. This implies that 
the explicit time dependencies in the radiation transfer 
and hydrodynamics (or in the given density profile and ve- 
locity field) can be neglected and that the time evolution 
of the nova atmosphere can be represented by a sequence 
of snapshot models. This assumption is valid since the 
hydrodynamic time scale is much greater than both the 
radiative time scale and also the excitation an d ionization 
equilibrium time scales l|Bath fc Shaviviri976|) . 

To account for LTE and NLTE line blanketing, the 
opacities of all important spectral lines have to be in- 
cluded. This also includes the handling of line blends, 
e.g., due to line broadening. Doppler broadening due to 
the large scale velocity field is handled through the co- 
moving frame radiation transport and the Lorentz trans- 
formation of the radiation field is used to obtain the ob- 
served spectrum in the Euler frame. PHOENIX includes a 
large number of NLTE lines as well as LTE background 
lines of species not treated in NLTE <|Hauschildtl 1993; 
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Hauschildt & Baronl Il995l, . It solves the line radiative 
transfer equations withou t using appro ximations such as 
the Sobolev method (see lGrinin1l200l|. and ref. therein). 
Therefore, all depth-dependent spectral line profiles have 
to be calculated in the Lagrangian frame. This in turn re- 
quires a large number of wavelength points (400,000 for 
the models presented in this paper). 

The generation of a model is very time consuming due 
to the complexity of the problem. The construction of 
detailed models is possible only throug h the use of par- 
allel algorithms in the PHOENIX code faauschild t et all 
Il997t iBaron k Hauschildt! Il99"3 lHauschildt et all 1200 1|) . 
Therefore, modern parallel computers are used to al- 
low more complex, more sophisticated, and more realistic 
models. 

2.1. Physical parameters of a nova model atmosphere 

There are several physical parameters of a nova model 
atmosphere. First there is the effective temperature T e g. 
Together with the reference radius r re f (radius at T s td = 1, 
where the optical depth grid is calculated at a standard 
wavelength, depending on the effective temperature of the 
model) it defines the bolometric luminosity of the nova 
via Lboi = 47rcrT^jr^ e f. This definition is only used as a 
convenient parameterization. 

There is ample evidence from modeling of the out- 
burst tb_Bt_the_jMya 1 _w^ by radiation pres- 
sure llHauschildt fc Starrfieldlfl994h . Therefore, the static 
Eddington luminosity L c( m cannot be used as an upper 
limit. Thus the luminosity of the nova has to be con- 
sidered as a free parameter or taken from hydrodynami- 
cal simulations. For the models presented here, we have 
assumed a bolometric luminosity of Lboi = 50, 000L©. 
This value implies a white dwarf mass close to the 
Chandrasekhar limit. In general, nova model atmospheres 
for earlier phases of the outburst have shown very little 
dependence of the spectrum on the bolometric luminosity 
(Hauschildt & Starrfieldlll995T) . This is also true for the 
models presented in this work and we parameterize our 
models by the effective temperature. 

The density profile of the nova atmosphere can be pa- 
rameterized as a power law of the form 



p{r) 



(1) 



The velocity profile follows from the requi rement that the 
mass loss rate M(r, t) is constant in space ( Bath fc Shavivl 
[[£3, i.e., 
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with two parameters N and v ou t- For novae, we use a typi- 
cal value of N — 3. v out is the outer velocity of the expand- 
ing shell, here we use the observed value of v out = 2500 km 
s _1 determined from the width of the line features of the 
X-ray spectra shown below. p out determines the outer 
edge of the envelope and is given by an outer pressure p out 



database 


X-ray lines 


all lines 


CHIANTI Ver. 3 


ps 12.650 


ps 45.000 


CHIANTI Ver. 4 


ps 22.350 


ps 74.500 


APED (total) 


ps 326.000 


w 463.200 


APED (distinct) 


ps 178.000 


w 310.400 



Table 1. Number of spectral lines from the atomic 
databases. Given are the number of X-ray lines (defined 
as lines with wavelengths A < 100 A) and the sum of all 
lines over the whole spectral range (up to A w 10 7 A). In 
the APED database entries for the same transition are 
possible. In the last row the number of unique lines are 
specified. 



assuming that the ideal gas equation is valid. Reasonable 
values lie in the range between p out = 10 -2 and 10 _1 dyn 
cm~ 2 . Together with N = 3 and r 2 in Eq. [3 this leads to 
a velocity field which is linear in r 

2.2. Additional physics included in PHOENIX for 
modeling novae in the SSS phase 

In order to model the X-ray spectra of novae with 
PHOENIX, we had to extend the opacity and line databases 
into the X-ray spectral range (from A ~ lA, equivalent 
to E f=s 12.4 keV). Several new atomic data and physi- 
cal processes have to be included in order to construct 
detailed nova atmosphere models with effective temper- 
atures in excess of 10 5 K. Previously the only database 
with X-ra y data avai l able i n PHO ENIX was CHIANTI 
Version 3 l)Dere et alJ Il997l l200l|) . For this work, we 
have extended PHOENIX to also include the CHIANTI 
Version 4 (CHIANTI4,|Yq 



mg et alJl2003l, and the APED 
(Astrophysical Plasma Emission Database, Version 1.3.1) 
1 databases. 

The new databases provide the following extensions 
compared to previous PHOENIX versions: 

— NLTE model atom data for highly ionized species or 
species with large ionization potential, including 

— many new spectral lines in the X-ray waveband down 
to 1 A 

— improved data for electron collision rates 

— new data for proton collision rates 

— better data for thermal bremsstrahlung 

All these improvements were used for the models pre- 
sented in this work. 

In Table ^ we give the number of atomic lines found 
in each of the databases. The APED database provides a 
good extension to the other databases. Note, that many 
of the lines in the databases are rather weak but, collec- 
tively, are important opacity sources. There are a number 
of transitions with multiple entries for the same line but 
different atomic data and it is not clear which entries are 
better than others. For this work, we simply use the first 



1 http://cxc.harvard.edu/atomdb/ 
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entry for each unique transition of the line databases and 
ignore all duplicate data. The effects of using different 
data where available will be investigated in a later paper. 

The values of the electron and proton collision rates 
for each ion depend significantly on the database. In ad- 
dition, the proton collision rates are only available in the 
CHIANTI4 and APED databases. All applicable collision 
rates are used in the rate equations of PHOENIX for ions 
whose transitions are treated in NLTE. A detailed inves- 
tigation of the sensitivity of the results on the collisional 
data will be presented in a subsequent paper. 

Thermal bremsstrahlung is calculated using the data 
for the ff-gaunt factors available in the databases. 
Previously, PHOENIX considered only contributions of H, 
Si, and Mg as bremsstrahlung processes (called 'old 
bremsstrahlung' in this paper). In the CHIANTI Version 
4 database, the ff-gaunt factors are available for elements 
up to Z = 30. These bremsstrahlung processes have been 
included in this work (called 'new bremsstrahlung' in this 
paper) which is important for highly ionized abundant el- 
ements. The nova model atmospheres are optically thick, 
therefore, free-free absorption as the inverse process to 
the bremsstrahlung emission is important and included in 
the models. Two ranges for the plasma temperature T and 
the ratio of photon to thermal energy u have to be dis- 
tinguished. For T > 10 6 K and u = > 1CT 4 (where v 
is the frequency of th e radiation ) the r elativistic ff-gaunt 
factors and the fits of lltoh et alJ l|2000j) are used. Outside 
this range we use the ff-gaunt factors o f lSutherlandl (l998) 
since they are valid over a large range of temperatures 
and frequencies. In both references, the ff-gaunt factors 
are given for elements up to Z = 28. However, it is pos- 
sible to enhance the fits up to Z = 30. These extensions 
are provided in the CHIANTI4 database and used in the 
model calculations presented here. 

2.3. Computation of hot nova models 

The models computed for this work are the first PHOENIX 
nova atmosphere models with effective temperatures 
T e s > 10 5 K. Due to the extreme conditions in nova atmo- 
spheres, it is better to start the model calculations from 
scratch rather than using scaled structures from our ear- 
lier nova modeling. The first step in this process is to 
calculate an LTE continuum model without any spectral 
lines for each temperature, typically starting from a grey 
temperature structure. In the next step, spectral lines are 
included in the LTE model to generate a model structure 
that includes LTE line blanketing. Next, NLTE models are 
constructed with more and more species treated in NLTE. 
It is important to treat the transitions of all important 
ionization stages of one element in NLTE simultaneously, 
because NLTE effects like over- and under-ionization af- 
fect the level populations of adjacent ionization stages. 
Therefore, large sets of NLTE species are needed in the 
model atmospheres. However, there are convergence prob- 
lems if the calculation includes too many 'fresh' NLTE 
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Fig. 1. The observed Chandra LETGS spectrum of nova 
V4743 Sgr. There are absorption lines and probably weak 
emission lines in the spectrum. The strongest lines are 
from the two hig hest ionisation stages of C, N, and O 
l|Ness et alJl2003|) . 

species at once. Therefore, we use a multi-stage process 
to reach the final models with the full set of NLTE species. 

Each step needs about 30 iterations and we calculate 
about 270 iterations in total for a model presented in this 
publication. The model includes about 7 x 10 3 b-f and 
about 100 x 10 3 b-b transitions in NLTE. The atmosphere 
is divided into 50 layers and distributed along an optical 
depth grid from T s u.min = 10~ 10 to T std _ max = 10 2 . 

3. X-ray observations of nova V4743 Sgr 

Nova V 4743 Sgr was first detected on 20 September 
2002 by lHaseda et all (2002) at my — 5. An accurate 
position is RA(J2000) = 19 ,1 01 m 09 s .38, DEC(J 2000) = 
-22°00'05".9 (determined bv iTanaka et alJl2002h . From 
the optical light curve the nova was classified as very 
fast. A velocity of 2400 km s _1 was estimated from the 
full width at half maximum (FWHM) of the Ha emis- 
sion line, ind icating that the e xpansion velocities are ~ 
1200 km s- MLvke et all ^2002) determined a distance of 
6.3 kpc based on infrared observations. The first X-ray 
observation of V4743 Sgr was obtained with the ACIS-S 
instrument on CHANDRA in 2002 November. It was faint 
at this time and did not appear to be in the SSS phase 
as yet. We then observed it for 24.7 ks with the HRC- 
S+LETG on 19 March 20 03 (OBSID 3775, iNess et alJ 
120031: IStarrfield et aljEool . 

The light curve for the first grating observation proved 
extremely interesting with both large amplitude oscilla- 
tions and, about two-thirds of the way through the ob- 
serv ation, a rapid de cline to almost no counts per sec- 
ond l|Ness et alJl200^1 . Two further HRC-S+LETGS ob- 
servations were carried out on 2003 July 18, and 2003 
September 26 (11.7ks: OBSID 3776 and 12.0 ks: OBSID 
4435). 

The spectrum from the first observation is shown in 
Fig.rj] The effective areas used to convert from ct s _1 to 
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flux are determined with the CIAO (Chandra Interactive 
Analysis of Observations) software package 2 , version 3.0. 
The spectrum is background subtracted. At wavelengths 
greater than ss 55 A, the spec trum is dominate d by second 
and higher dispersion orders l|Ness et al.ll2003[) . and these 
wavelengths will not be considered in our analysis. 

An examination of the spectrum shows that it is not 
a black-body but resembles a stellar atmosphere with 
deep absorption features and, possibly, some weak emis- 
sion lines. This spe ctrum onl y slightly resembles that 
of either CAL 83 JPaerels et alJ Ejblj) or V1494 Aql 
llStarrfield et alJl200lh . The strongest lines ar e from the 
two h ighest ionisation stages of C, N, and O llNess et, alJ 
l2003h . 

An extensive analysi s of the ob s ervati on from March 
has been carried out bv lNess et alJ l|2003j) . A remarkable 
feature in the light curve from that observation is a strong 
de cline of the X-ray emission between 13 and 20 ks. Fig. 1 ^ 
of iNess et al.l 1)200311 shows the light curve taken during 7 
this observation. It is oscillating at high amplitude with „ 
a period of ~ 22 min. Moreover, the oscillations exhibit E 
spectral variability and a correlated change in hardness 
ratio, defined by [h — s)/(h + s), with h being the flux in ^ 
19 A < A < 30 A and s being the flux in 30 A < A < 50 * 
A. There is a strong period present in the data, with at 
least two harmonic overtones, which reflect the complex 
nature of the variations. This is in line with the X-ray os- 
cillations of V1494 Aql which w ere also mult i perio dic and 
were interpreted as pulsation bv lDrake et al.ll|2003h . If the 
oscillations in V4743 Sgr (which exhibited a much larger 
amplitude than those in V1494 Aql) are interpreted as a 
rotational modulation then they imply a major asymme- 
try in the X-ray emitting region. 

The spect r um of the last 5.4 ks can be seen in Fig. 4 
of INess et al.1 l(2003j) . There are strong emission lines but 
there is no continuum present so we cannot discern if ab- 
sorption lines are present. In any case, the peak count rate 
at this time is far lower than in any of the X-ray obser- 
vations where a strong continuum is present. A possible 
scenario for the decline at 13 ks involves an eclipse of the 
X-ray emitting WD. But this does not explain either the 
long duration of the decline or the softening during the 
decline. The emission spectrum at minimum is probably 
produced by the ejected shell and not a hot stellar pho- 
tosphere. Therefore, we do not consider it any further in 
this paper. 

4. Results 

The PHOENIX code package was modified to calculate 
model atmospheres and synthetic spectra for the late 
phases of classical novae where most of the energy is 
emitted in X-rays. Our models have effective tempera- 
tures of T e ff = 10 5 — 10 6 K and a bolometric luminosity of 
Lbo\ = 50, OOOLq. In the generation of the first test mod- 
els the 'old bremsstrahlung' was used. A typical NLTE 

2 |http : // exc . h arvard . edu/ ciao/ 
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Fig. 2. Upper: Model spectrum calculated with the 'old 
bremsstrahlung' continuum. Because these are test mod- 
els, there are only spectral lines from H, He, C, and 
N. Lower: Model spectrum calculated with the 'new 
bremsstrahlung' continuum and spectral lines of H, He, 
C, N, O, and Fe. There are mainly absorption lines in the 
spectrum. 

model spectrum can be seen in the upper panel of Fig. 
There are only a few weak absorption lines and the other 
spectral lines are in emission. Clearly, the observed spec- 
tra (Fig. ^) are not well represented by models using the 
'old bremsstrahlung'. 

In calculating the other grids of models the 'new 
bremsstrahlung' was used (Fig. [21 lower panel), where 
there are mostly absorption lines in the spectrum. 
Considering the interstellar and circumstellar absorption 
by hydrogen, the shape of the continuum model spectrum 
is close to what is observed (Fig. [3J| . The agreement be- 
tween the continuum of the model and the observation is 
acceptable, whereas there are differences in the spectral 
lines. 

The best fit to the spectrum of nova V4743 is found 
at a temperature of T e ff = 5.8 x 10 5 K. To get the cor- 
rect slope for the pseudo-continuum (continuum formed 
by numerous overlapping lines), a suitable value of the 
hydrogen column density has to be used. A value of 
rih — 4.0 x 10 21 cm~ 2 for the hydrogen column den- 
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T eff = 5.8 10 5 K, n h = 4.0 10 21 cm" 
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Fig. 3. NLTE Model spectrum (calculated with the 'new 
bremsstrahlung') with solar abundances plotted over the 
first 15 ks of the March observation. There is a very 
good agreement. A best-fit is provided with an effective 
temperature of T off = 5.8 x 10 5 K (L hol = 50,000L o , 
v out = 2500 km/s, p xet = 10- 2 dyn/cm 2 , N = 3). The 
absorption lines and the continuum absorption are too 
weak in the model. Interstellar absorption was taken into 
account (n h = 4.0 x 10 21 cm~ 2 ). 



sity provides the best fit, in contrast to the value o f 
rih = 1.41 x 10 21 cm -2 given bv lDickev fc Lockmanl l)l990|) 
for the ISM. This discrepancy can be explained by assum- 
ing additional absorption from the circumstellar shell, not 
included in the value from the literature. 

Close inspection of the measured spectrum reveals that 
some spectral lines are not reproduced well or are miss- 
ing in the model spectrum. This is because we have used 
only solar abundances in the model and have not en- 
hanced abundances of, for example the CNO elements, 
as is generally observed in novae and predicted by theory. 
In the best model spectrum with solar abundances, all 
absorption lines are too weak and there is too much emis- 
sion around 24 A (ss 517 eV). Increasing the abundances 
should increase the absorption. We will do this in a suc- 
ceeding publication. There is an oxygen edge around 24 A. 
Increasing the O-abundance should add more opacity and 
could provide a better fit. The computation of models with 
abundances deviating from solar needs much more com- 
putational work and has not yet been done. With such 
models it will be possible to estimate the chemical abun- 
dances of the nova ejecta. 

The lack of lines in the computed spectrum may also 
be caused by using only line data from ions of H, He, C, 
N, and O and from Fe XXI-XXVI in the model calcula- 
tions. If the concentration of the partial pressures of an 
ion is sufficiently large, all available transitions between 
all available levels of these ions are treated in NLTE. This 
procedure selects the last three ionization stages of C, N, 
and O to be included in the NLTE calculation and the 
others of these elements are treated in LTE. 



Departure coefficients of H 
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T std 

Departure coefficients of C IV — V 




Departure coefficients of N V-V 
10 6^ 
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Departure coefficients of VI — V 
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Fig. 4. Departure coefficients bi from the best-fit model 
for the first 5 levels of the H-, C-, N-, and O-ions consid- 
ered in these models. There are departure coefficients with 
values of 6 magnitudes larger or smaller than 1. For all el- 
ements the departures from LTE are strong. As expected, 
the departures are strongest for outer layers and there is 
no departure from LTE at large optical depths (pi = 1). 
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Fig. 5. Upper: Radius of the nova atmosphere plotted 
against the standard optical depth. The nova atmosphere 
is very extended, with a radius of r ~ 10 12 cm at r s td = 

io- 12 . 

Lower: Temperature distribution in the atmosphere plot- 
ted against the standard optical depth. There is a wide 
variety of temperatures in the atmosphere. In the inner 
layers the temperature is several 10 6 K high. 



For each ion it can be decided which atomic database 
will be used. Because many X-ray spectral lines are 
required, the CHIANTI4 or APED databases are used 
(Tabic [TJ. From a rough comparison, the spectra calcu- 
lated with each database are only slightly different. The 
shape of the continuum and the number of strong lines are 
the same. The quality of the databases has to be deter- 
mined for an accurate abundance analysis. For the models 
presented in this paper we have used the APED database. 

The main results from the models are that the X-ray 
emission is dominated by thermal bremsstrahlung from 
the atmosphere surrounding the WD, and that the hard 
spectral range of A < 29 A (E > 428 eV) is dominated 
by iron and nitrogen absorption. As expected, the depar- 
tures from LTE are large in the nova ejecta. To examine 
the strength of NLTE effects, the departure coefficients 
can be plotted against the standard optical depth r st( i- 
The departure coefficient hi is the ratio of the occupation 
numbers in level i calculated in NLTE to calculations in 
LTE for a specific ion. The stronger the departure from 



Partial pressures of FelX — XXV 




10 u 10 10 10 10 10 10 



Fig. 6. Partial pressures of FelX-XXVII. The highest ion- 
ization stages of iron exist at the deepest layers of the nova 
atmosphere. In the outer layers FeXVII is most abundant. 
The partial pressures of Fel to FeVIII are very small and 
lie below the range of the y-axis. Only ions from FeXV are 
abundant enough in the line forming region and therefore 
important. 

LTE the more different the value of the departure coeffi- 
cient from 1. In the outer layers of nova atmospheres the 
departure coefficients can differ up to six dex from 1 in 
both directions. This can be seen in Fig. 0] for different 
elements. Spectral lines are formed in regions where the 
departure coefficients are very different from 1. Therefore 
it is very important to treat all atomic transitions in full 
NLTE. 

In addition to strong departures from LTE, there are 
extreme physical conditions in the nova ejecta. The atmo- 
sphere is very extended with a radius of r « 10 12 cm at 
T s td = 10~ 12 (Fig.0 upper panel). There is a wide range 
of electron temperatures (up to several 10 6 K, Fig.0 lower 
panel) and hence several of the highest ionization stages 
of an element are simultaneously present in the ejecta. For 
example the partial pressures of iron ions (Fe IX-XXVII) 
are plotted against the standard optical depth in Fig. El A 
similar effect has alr eady been observed in the IUE anal- 
ysis of nova spectra l)Hauschildt et al.lll992|) . 

5. Conclusion and outlook 

We have modified the PHOENIX code to calculate model 
atmospheres and synthetic spectra for the late X-ray 
phase of classical novae. The agreement between the 
model and the observed spectrum of nova V4743 Sgr is 
good in the continuum for an effective temperature of 
T e g = 5.8 x 10 5 K if the 'new bremsstrahlung' (see sec- 
tion 122} is used which is physically more reasonable. So 
far, we only have computed model atmospheres with so- 
lar abundances for computing time reasons and the fit to 
the spectral lines needs to be improved by models with 
non-solar abundances. The hydrogen column density of 
rih = 4.0 x 10 21 cm~ 2 providing the best fit is high com- 
pared to — 1.41 x 10 21 cm~ 2 from lDickev fe Lockmanl 
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(1990), because of additional absorption by the circum- 
stellar shell. 

X-ray continuum emission is dominated by thermal 
bremsstrahlung and its inverse processes, while bound- 
free emission is much weaker in the nova atmosphere. In 
the hard spectral range of A < 29A (E > 428eV) iron 
and nitrogen absorption are dominant. As expected, nova 
ejecta show strong departures from LTE. They are very 
extended and several of the highest ionization stages of all 
elements are simultaneously present. 

In future work additional physical processes absorbing 
and emitting X-rays will be implemented in the code. For 
example there are data for dielectronic recombination and 
autoionization, and there are better data for free-bound 
coefficients in the APED and CHIANTI4 databases. 

Models with non-solar abundances will be calculated 
to obtain a better fit for the line strengths and the con- 
tinuum absorption and to analyze the abundances of the 
nova ejecta. To determine reliable chemical abundances 
more elements have to be included in the NLTE calcula- 
tions. 

The data quality of the atomic databases as well as 
differences between the CHIANTI4 and APED database 
will be checked, e. g. by examining the line profiles. By 
comparing model and observation, it has to be determined 
which data of the individual ions reproduce the spectral 
lines as accurately as possible. This is necessary for an 
accurate abundance analysis. 

The evolution of the X-rays should be examined, by 
modeling the spectra of other observations of novae with 
the CHANDRA satellite. It is important to test if it is 
possible to model the spectrum of novae in the emission 
line phase, when nuclear burning has switched off and the 
exp anding shell is optica lly thin, like the spectrum of V382 
Vel (Bur witz et al . 2002). For this work more observations 
with CHANDRA are necessary. 
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